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Abstract  Mathematical  modeling  of  an  Electric-Swing  Ad¬ 
sorption  (ESA)  system  (adsorption  cycle  with  electrother¬ 
mal  desorption  step,  performed  by  direct  heating  of  the  ad¬ 
sorbent  particles  by  passing  electric  current  through  them), 
with  annular,  radial-flow,  cartridge-type  fixed-bed  and  in¬ 
vessel  condensation,  is  performed  by  using  Comsol  Multi¬ 
physics™  software.  Three  multiphysics  models  are  built,  in 
order  to  describe  three  stages  of  a  compete  ESA  cycle:  ad¬ 
sorption,  electrothermal  desorption  before  the  start  of  con¬ 
densation  and  electrothermal  desorption  with  in- vessel  con¬ 
densation.  In  order  to  describe  the  complete  ESA  cycle  the 
models  for  the  three  stages  are  integrated,  by  using  a  combi¬ 
nation  of  Comsol  Multiphysics™  and  Matlab™.  The  mod¬ 
els  were  successfully  used  for  simulation  of  separate  stages 
of  the  process  and  of  the  complete  ESA  cycles,  as  well  as 
for  investigation  of  the  influences  of  the  main  operational 
parameters  on  the  process  performance. 
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Abbreviations 

a  (m2/m3) 

Specific  surface  area 

b  (K-1) 

Temperature  coefficient  of  the  bed 
electrical  resistivity 

C  (mol/mol) 

Adsorbate  concentration  in  the 
gas  phase 

C*(mol/mol) 

Adsorbate  concentration  in  the 
gas  phase  in  equilibrium  with  the 
solid  phase 

Cbreak  (mol/mol) 

Breakthrough  concentration 

Csai  (mol/mol) 

Saturation  concentration 

Cpg  (J/mol/K) 

Specific  heat  capacity  of  the  inert 
gas 

Cpi  (J/mol/K) 

Specific  heat  capacity  of  liquid 
adsorbate 

cps  (J/g/K) 

Heat  capacity  of  the  solid  phase 

cpv  (J/mol/K) 

Specific  heat  capacity  of  gaseous 
adsorbate 

Dm  (mol/cm/s) 

Mass  transfer  dispersion 
coefficient 

d’!r  (W/K/cm) 

Heat  diffusivity  of  the  gas  phase 

D\,s  (W/K/cm) 

Heat  diffusivity  of  the  solid  phase 

E  (J/mol) 

Adsorption  energy  of  the 
adsorbate  (D-R  equation) 

8  (cm/s2) 

Gravitation  constant 

G  (mol/s) 

Flow  rate  of  the  inert  gas 

H  (cm) 

Bed  axial  dimension  (Fig.  2) 

h  i  (cm) 

Adsorber  dimension  (Fig.  2) 

hb  (J/cm2/K) 

Solid  to  gas  heat  transfer 
coefficient  within  the  bed 
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hsi  (J/cm2/K) 

Heat  transfer  coefficient  from  the 
solid  phase  to  the  gas  phase  in  the 
central  tube(s) 

hs 2  (J/cm2/K) 

Heat  transfer  coefficient  from  the 
solid  phase  to  the  gas  phase  in  the 
annular  tube 

hwg  (J/cm2/K) 

Gas  to  ambient  heat  transfer 
coefficient  (heat  losses) 

J  (A/cm2) 

Current  density 

Jcond  (mol/cm2/s) 

Condensation  flux 

k  (cm2) 

Asorbent  bed  permeability 

km  (mol/cm2/s) 

Mass  transfer  coefficient  in  the 

adsorbent  bed 

km  i  (mol/cm2/s) 

Mass  transfer  coefficient  between 
the  adsorbent  bed  and  the  gas  in 
the  central  tube 

km 2  (mol/cm2/s) 

Mass  transfer  coefficient  between 
the  adsorbent  bed  and  the  gas  in 
the  annular  tube 

Lcond  (mol/s) 

Flow-rate  of  the  condensed  liquid 

Lcond  (mol) 

Total  amount  of  the  condensed 
liquid 

P  (Pa) 

Gas  pressure 

Pa  (Pa) 

Ambient  pressure 

Pc  (Pa) 

Critical  pressure 

p°  (Pa) 

Adsorbate  saturation  pressure 

Qel  (W) 

Rate  of  heat  generation  (equal  to 
electric  power) 

Qel  (J) 

Heat  generation  (equal  to  electric 
energy) 

q  (mol/ g) 

Adsorbate  concentration  in  the 
solid  phase 

Rg  (J/mol/K) 

Gas  constant 

r  (cm) 

Radial  coordinate 

r\  (cm) 

Adsorber  dimension — radius  of 
the  central  tube  (Fig.  2) 

r2  (cm) 

Adsorber  dimension  (Fig.  2) 

r3  (cm) 

Adsorber  dimension  (Fig.  2) 

r4  (cm) 

Adsorber  dimension — radius  of 
the  adsorber  vessel  (Fig.  2) 

Sw  (cm2) 

Surface  are  of  the  outer  wall  of 

the  adsorber  vessel 

Ta  (K) 

Ambient  temperature 

Tg  (K) 

Gas  phase  temperature 

Tc  (K) 

Critical  temperature 

Tr  (K) 

Referent  temperature 

Ts  (K) 

Solid  phase  temperature 

Tsw  (K) 

Switch  temperature  (from 
desorption  to  adsorption) 

Tw  (K) 

Wall  temperature 

t  (s) 

Time 

U(V) 

Electric  potential 

u()  (V) 

Supply  voltage 

u  (cm/s) 

Radial  superficial  gas  velocity 

VPA ,  VPB ,  VPC ,  VPD 

Wagner  constants  (Wagner 

v  (cm/s) 

equation) 

Axial  superficial  gas  velocity 

Wo  (cm3/g) 

Total  volume  of  micropores  (D-R 

z  (cm) 

equation) 

Axial  coordinate 

Greek  letters: 

ap 

Purification  factor 

Us 

Separation  factor 

A  Hods  (J/mol) 

Heat  of  adsorption 

A  Hcond  (J/mol) 

Heat  of  condensation 

£b 

Bed  porosity 

M  (Pa/s) 

Dynamic  viscosity 

p  (£2cm) 

Electric  resistivity 

po  (£2cm) 

Electric  resistivity  at  referent 

pg  (mol/cm3) 

temperature  Tr 

Gas  phase  density 

Pb  (g/cm3) 

Adsorbent  bed  density 

pA  (g/cm3) 

Adsorbate  density 

X  (s) 

Time  period 

Subscripts: 

A 

Adsorption 

b 

Bed 

D 

Desorption 

8 

Gas 

in 

Inlet 

out 

Outlet 

P 

Previous  (initial) 

r 

In  the  radial  (r)  direction 

s 

Solid  phase 

ct 

Central  tube 

as 

Annular  space 

z 

In  the  axial  (z)  direction 

Special  characters: 

0 

Time  average  value 

1  Introduction 

Electric  (or  electrothermal)  swing  adsorption  (ESA)  is  a  rel¬ 
atively  new  name  (Yu  et  al.  2007;  Sullivan  et  al.  2004a)  for 
an  essentially  temperature- swing  adsorption  (TSA)  cycle  in 
which  the  adsorbent  material  is  heated  by  passing  electric 
current  through  it  (owing  to  Joule  effect). 

The  idea  about  regeneration  of  used  adsorbent  beds  by 
direct  heating  by  electric  current  was  first  published  in  the 
1970s  (Fabuss  and  Dubois  1970).  The  desorption  process 
based  on  this  principle  was  later  named  electrothermal  des¬ 
orption  (Petkovska  et  al.  1991)  and  it  was  recognized  as  a 
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prospective  way  to  perform  desorption  steps  of  TSA  cycles 
(Petkovska  et  al.  1991). 

The  main  differences  between  electrothermal  desorption 
and  classical  techniques  of  thermal  regeneration  of  used  ad- 
sorpbents  (e.g.  heating  by  hot  gas  or  steam)  are  (Sullivan 
2003;  Petkovska  and  Mitrovic  1994a,  1994b): 

•  Energy  efficiency  is  higher  than  when  heating  by  steam 
or  inert  gas,  because  the  energy  is  delivered  directly  to 
the  adsorbent,  thus  minimizing  the  energy  expended  for 
heating  of  the  vessel  and  ancillary  equipment. 

•  Heating  rate  of  the  adsorbent  is  not  limited  by  heat/mass 
transfer  rates  between  the  carrier  gas  and  the  adsor¬ 
bent,  because  the  energy  is  directly  applied  to  the  ad¬ 
sorbent  rather  than  being  delivered  by  the  carrier  gas 
stream. 

•  Effluent  concentration  of  adsorbate  can  be  maximized, 
because  the  purge  gas  flow  rate  is  controlled  indepen¬ 
dently  from  the  power  applied  to  the  adsorbent. 

•  Unlike  steam  regeneration,  no  water  is  introduced  into 
the  system,  which  is  difficult  and  costly  to  separate  from 
water-miscible  adsorbates,  and  which  can  cause  corro¬ 
sion. 

•  The  directions  of  heat  and  mass  transfer  are  identical  in 
the  case  of  electrothermal  desorption,  and  opposite  in  the 
case  of  conventional  thermal  desorption  by  hot  gas  or 
steam.  This  can  influence  the  process  kinetics  and  dy¬ 
namics. 

It  was  shown  that  electrothermal  desorption  has  some  advan¬ 
tages  over  conventional  methods,  regarding  adsorption  ki¬ 
netics  and  dynamics  (Petkovska  and  Mitrovic  1994a,  1994b) 
and  energy  efficiency  (Petkovska  et  al.  1991;  Sullivan  2003). 

Most  applications  of  electrothermal  desorption  have  been 
based  on  fibrous  activated  carbon  materials  (Petkovska  et  al. 
1991;  Rood  et  al.  2002;  Sullivan  et  al.  2001,  2004a,  2004b; 
Subrenat  and  Le  Cloirec  2004,  2006;  Moon  and  Shim  2006; 
Luo  et  al.  2006),  although  granular  activated  carbons  (Sny¬ 
der  and  Leesch  2001;  Sushchev  et  al.  2002;  Luo  et  al.  2006) 
and,  recently,  carbon  monoliths  (Luo  et  al.  2006;  Yu  et  al. 
2002,  2007),  have  also  been  used. 

Up  to  now,  different  applications  (mostly  on  laboratory 
scale)  of  regeneration  of  used  adsorbents  by  electrother¬ 
mal  desorption  have  been  reported:  air  purification,  with 
emphases  on  removal  of  hazardous  volatile  organic  com¬ 
pounds  (Baudu  et  al.  1992;  Azou  et  al.  1993;  Sullivan  et 
al.  2001,  2003;  Subrenat  and  Le  Cloirec  2006),  gas  separa¬ 
tion  (Moon  and  Shim  2006)  and  solvent  recovery  (Bathen 
et  al.  1997;  Lordgooei  et  al.  1996;  Snyder  and  Leesch  2001; 
Sullivan  2003).  In  the  last  couple  of  years,  some  applica¬ 
tions  of  electrothermal  desorption  on  industrial  scale  have 
also  been  reported  (Bathen  2002;  Subrenat  and  Le  Cloirec 
2004,  2006). 

The  subject  of  this  investigation  is  mathematical  mod¬ 
eling  of  an  ESA  system  based  on  assembly  of  annular, 


radial-flow,  cartridge-type  fixed-beds,  with  in- vessel  con¬ 
densation  (Rood  et  al.  2002;  Sullivan  2003;  Sullivan  et  al. 
2004a),  which  has  been  developed  for  removal  and  recov¬ 
ery  of  hazardous  volatile  organic  compounds.  This  system 
involves  a  large  number  of  coupled  processes  and  phenom¬ 
ena:  fluid  flows  through  tubes  and  porous  beds,  adsorption, 
mass  transfer,  heat  transfer,  heat  generation  associated  with 
electric  current  flow,  and  condensation,  all  in  a  rather  com¬ 
plex  geometry.  Such  a  system  is  very  demanding  from  the 
modeling  point  of  view  and  represents  a  good  candidate  for 
multiphysic  modeling.  This  manuscript  presents  the  proce¬ 
dure  and  results  of  mathematical  modeling  of  this  system  by 
using  Comsol  Multiphysics™,  a  software  tool  specialized 
for  multiphysics  modeling. 

2  Description  of  the  modeled  system 

The  adsorber  of  the  ESA  system  treated  in  this  work  con¬ 
tains  one  cartridge-type  adsorption  bed.  It  is  schematically 
shown  in  Lig.  la.  The  cartridge,  which  is  schematically 
shown  in  Lig.  lb,  is  formed  as  a  cylindrical  roll  of  activated 
carbon  fiber  cloth,  spirally  coiled  around  a  porous  central 
pipe  used  for  introducing  the  gas  stream  (the  stream  that  has 
to  be  purified,  during  adsorption,  or  an  inert  gas  stream,  dur¬ 
ing  desorption).  The  gas  flow  through  the  adsorption  bed  is 
in  the  radial  direction. 

During  the  desorption  step,  a  low  flow-rate  stream  of 
pure  inert  is  introduced  into  the  adsorber  and,  at  the  same 
time,  electric  current  is  passed  through  the  adsorbent  bed 
in  the  axial  direction,  causing  heat  generation,  heating  of 
the  bed  and  desorption.  The  desorbed  vapor  is  carried  by 
the  radially  flowing  gas  stream  towards  the  adsorber  wall. 
When  the  vapor  concentration  near  the  wall  becomes  equal 
to  the  saturation  concentration  corresponding  to  the  wall 
temperature,  the  vapor  is  condensed  and  drained  out  of 
the  system  through  the  funnel  shaped  bottom  outlet  (see 
Lig.  la). 

The  switching  between  adsorption  and  desorption  is  per¬ 
formed  when  the  outlet  concentration  reaches  some  prede¬ 
fined  breakthrough  value.  In  order  to  avoid  thermal  degra¬ 
dation  of  the  adsorbate,  the  desorption  step  is  ended  by 
turning  the  electric  heating  off  when  the  temperature  in  the 
adsorbent  bed  reaches  some  predefined  maximal  tempera¬ 
ture. 

The  following  assumptions  were  used  in  setting  up  the 
models: 

•  The  adsorbent  bed  is  treated  as  homogeneous,  with  uni¬ 
form  adsorbent  porosity  and  density. 

•  The  mass  and  heat  transfer  resistances  on  the  particle 

scale  are  neglected. 

•  The  fluid  phase  is  treated  as  an  ideal  gas  mixture  of  an 

inert  and  a  single  adsorbate. 
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GAS AND 

CONDENSATE  OUT 

(a) 


Fig.  1  Schematic  representation  of:  (a)  the  adsorber  with  car¬ 
tridge-type  adsorbent  bed  and  in-vessel  condensation;  and  (b)  the  an¬ 
nular,  radial-flow  cartridge 


•  The  electric  resistivity  of  the  adsorbent  is  temperature 
dependent.  Linear  temperature  dependence  is  assumed, 
based  on  experimental  results  reported  in  (Sullivan  2003). 

•  All  other  physical  parameters  and  coefficients  are  consid¬ 
ered  as  constants. 

•  The  condensation  at  the  adsorber  wall  is  drop  wise.  This 
assumption  is  based  on  experimental  observations  re¬ 
ported  in  Sullivan  (2003). 

•  The  volume  of  the  condensed  liquid  is  neglected,  i.e.,  it  is 
assumed  that  the  liquid  drops  don’t  influence  the  gas  flow. 

•  The  heat  resistance  and  heat  capacity  of  the  adsorber  wall 
are  neglected,  so  the  wall  temperature  is  regarded  as  equal 
to  the  temperature  of  the  environment. 

•  Heat  transfer  by  radiation  is  neglected. 

•  The  electric  power  during  electrothermal  desorption  is 
supplied  under  constant  voltage  conditions. 


3  Multiphysics  models  for  separate  steps  of  the  ESA 
cycle 

Analysis  of  the  adsorber  defined  in  the  previous  section 
shows  that  it  consists  of  three  parts:  the  central  tube  used  for 
introducing  the  gas,  the  adsorbent  bed  and  the  annular  space 
between  the  adsorbent  bed  and  the  outer  adsorber  wall.  Also, 
three  consecutive  steps  of  the  complete  ESA  cycle  with  in¬ 
vessel  condensation  can  be  defined: 

•  Adsorption; 

•  Desorption  without  condensation; 

•  Desorption  accompanied  with  condensation. 

The  last  two  steps  are  parts  of  the  desorption  half-cycle. 

The  mathematical  model  of  the  system  consists  of  the 
momentum,  mass  and  heat  balances  for  different  parts  of  the 
system  (the  central  tube  of  the  cartridge,  the  adsorbent  bed 
and  the  annular  space  around  the  cartridge)  and  for  the  gas 
and  solid  phases,  as  well  as  the  balance  equation  defining 
electric  current  flow  and  heat  generation  owing  to  Joule’s 
effect. 

In  addition,  the  Dubinin-Radushkevich  (D-R)  equation 
(Dubinin  1989)  is  used  to  calculate  the  adsorption  equilib¬ 
rium  in  this  system. 

In  principle  all  three  steps  of  the  ESA  cycle  can  be  de¬ 
scribed  by  the  same  set  of  model  equations.  Nevertheless, 
the  heat  balance  for  the  solid  phase  and  the  boundary  condi¬ 
tions  corresponding  to  different  steps  can  be  different,  so  it 
was  necessary  to  build  three  different  models.  Let  us  name 
the  models  for  adsorption,  desorption  and  desorption  with 
condensation:  Model_A,  Model_D  and  Model_DC,  respec¬ 
tively. 

In  order  to  reduce  the  model  complexity,  and  using  the 
fact  that  the  adsorber  is  practically  axially  symmetrical,  the 
Comsol  Multiphysics  models  were  built  in  the  axially  sym¬ 
metrical  2D  space.  The  model  geometry,  which  is  identical 
for  all  three  models,  is  shown  in  Fig.  2.  As  can  be  seen  from 
this  figure,  the  model  is  defined  for  only  one  half  of  the 
adsorber,  and  the  other  half  is  axially  symmetrical.  In  that 
way,  the  number  of  the  needed  mesh  elements  is  reduced  to 
half,  which  dramatically  reduces  the  computational  time.  In 
Fig.  2,  we  also  define  the  main  adsorber  dimensions  and  the 
sub-domains.  Based  on  the  structure  of  the  adsorber,  three 
sub-domains  are  defined  (Fig.  2): 

•  SI — the  central  tube; 

•  S2 — the  adsorbent  bed; 

•  S3 — the  annular  space  between  the  cartridge  and  the  outer 
adsorber  wall. 

In  the  Comsol  Multiphysics  working  environment,  the 
model  equations  can  be  defined  by  choosing  the  appropriate 
application  modes  from  a  menu.  The  following  application 
modes  were  used  in  our  models: 
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Fig.  2  Definition  of  the  model 
geometry  and  sub-domains 
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no  condensation  at  the  adsorber  wall  is  present,  zero  mass 
flux  at  the  wall  is  assumed,  while  the  heat  flux  is  defined 
by  heat  losses  (35).  On  the  other  hand,  in  the  case  of  con¬ 
densation,  the  boundary  condition  for  the  mass  balance 
is  defined  by  constant  concentration  (equal  to  the  satu¬ 
ration  concentration  corresponding  to  the  wall  tempera¬ 
ture),  while  an  additional  term  defining  the  heat  of  con¬ 
densation  is  added  to  the  equation  defining  the  boundary 
condition  for  the  heat  balance  (36). 

The  seven  application  modes  listed  here  define  the  ve¬ 
locities,  concentrations  and  temperatures  in  all  three  sub- 
domains  and  both  phases,  as  well  as  the  electric  potential 
and  resistance  of  the  adsorbent  bed,  as  functions  of  time  and 
space  coordinates.  Nevertheless,  it  is  recommendable  to  de¬ 
fine  and  calculate  two  additional  parameters  which  charac¬ 
terize  the  ESA  process:  the  mass  of  the  condensed  liquid  per 
cycle  and  the  used  electric  energy  per  cycle. 


1.  N on-isothermal  flow — for  defining  momentum  balances 
and  fluid  flows  in  tubes  (in  our  models  defined  in  sub- 
domains  SI  and  S3); 

2.  Brinkman  Equations — for  defining  momentum  balances 
and  fluid  flows  in  packed  beds  (in  sub-domain  S2); 

3.  Convection  and  Diffusion — for  defining  mass  transfer  in 
flowing  fluids  (in  all  three  sub-domains); 

4.  Diffusion — for  defining  mass  transfer  in  stagnant  phases 
(in  sub-domain  S2); 

5.  Convection  and  Conduction — for  defining  heat  transfer 
in  flowing  fluids  (in  all  three  sub-domains); 

6.  Heat  transfer  by  Conduction — for  defining  heat  transfer 
in  stagnant  phases  (in  sub-domain  S2); 

7 .  Conductive  Media  DC — for  defining  electric  current  flow 
and  heat  generation  by  Joule’s  effect  (in  sub-domain  S2). 

In  principle,  the  Conductive  Media  DC  application  mode 
needs  to  be  defined  only  in  Model_D  and  Model_DC,  corre¬ 
sponding  to  the  parts  of  the  ESA  cycle  when  the  electricity  is 
on.  Nevertheless,  in  order  to  integrate  the  three  models  into 
a  model  of  the  complete  ESA  cycle,  this  application  mode 
was  used  in  all  three  models. 

The  complete  mathematical  model  of  our  system  is  too 
cumbersome  for  presenting  it  here.  The  model  equations  are 
listed  in  Appendix  1 .  Here  we  only  point  out  the  main  dif¬ 
ferences  between  the  three  models: 


3.1  Calculation  of  the  mass  of  the  condensed  liquid 

The  flux  of  the  condensed  liquid  is  calculated  from  the  gas 
concentration  gradient  at  the  adsorber  wall: 

(1) 

r=r4 

This  flux  is  a  function  of  the  axial  position  z  and  time  t  (r 
is  the  radial  coordinate,  r\  the  radius  of  the  adsorber  vessel 
(see  Fig.  2),  C  the  gas  concentration,  and  Dm  r  the  radial  dif¬ 
fusion  coefficient  of  the  gas  phase).  The  condensation  rate  is 
calculated  by  integration  of  this  flux  over  the  whole  surface 
of  the  adsorber  wall,  and  is  a  function  of  time  only: 

r  Sw 

Lcond(t)  =  /  Jcond(Zi  t)dS .  (2) 

JO 

The  total  amount  of  the  condensed  liquid  is  obtained  by  in¬ 
tegrating  the  condensation  rate  over  time: 

r ? cond  . 

Lcond  —  I  Lcond(t)dt.  (3) 

Jo 

Sw  is  the  total  surface  area  of  the  outer  adsorber  wall  and 
rcond  the  total  time  of  condensation. 

3.2  Calculation  of  the  used  electric  energy 


nd(Zi  0  —  Dm?r  „ 


dC 

dr 


1 .  The  energy  balance  for  the  solid  phase  ((20)  from  Appen¬ 
dix  1)  contains  a  term  defining  the  Joule  heat  generation 
(8  Q ei /dV).  For  Model_A  this  term  is  equal  to  zero. 

2.  The  boundary  conditions  for  the  mass  and  heat  bal¬ 
ances  corresponding  to  the  adsorber  wall  are  different  in 
Model_A  and  Model_D  on  one  hand  (the  parts  of  the 
process  without  condensation)  and  Model_DC,  on  the 
other  hand  (the  part  of  the  process  with  condensation).  If 


The  used  electric  energy  is  equal  to  the  Joule  heat  produced 
in  the  adsorbent  material.  This  energy  per  unit  volume  and 
unit  time  is: 


SQel  _  1  ((d_y_f  (duf\ 
dV  p  V  \  9r  /  \  dz  )  ) 


(4) 


where  p  is  the  electric  resistivity,  which  is  temperature  de¬ 
pendent,  and  U  the  electric  potential. 
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The  used  electric  power  is  obtained  by  integration  over 
the  adsorbent  bed  volume: 


QelO )  — 


and  the  total  used  energy  during  desorption,  by  integration 
over  time: 

P  Tdes 

Qel  —  /  Qel(t)dt.  (6) 

Jo 

Tdes  is  the  total  desorption  time  and  //,  r 2  and  r 3  the  adsorber 
dimensions  defined  in  Fig.  2. 


4  Model  of  the  complete  ESA  cycle 

In  order  to  simulate  the  complete  ESA  cycle,  the  three 
Comsol  Multiphysics  models  defined  in  the  previous  sec¬ 
tion  have  to  be  run  in  a  sequence,  with  automatic  switch¬ 
ing  from  one  model  to  the  next  one  when  certain  predefined 
conditions  are  met.  The  following  switches  need  to  be  per¬ 
formed: 

•  Switching  from  Model_A  to  Model_D.  This  switch  corre¬ 
sponds  to  the  end  of  adsorption  and  start  of  a  new  desorp¬ 
tion  step.  It  is  performed  automatically  when  the  outlet 
concentration  reaches  a  certain  predefined  breakthrough 
value  C break-  The  outlet  concentration  is  obtained  as  the 
average  value  across  the  outlet,  although  the  simulation 
results  showed  that  the  gas  concentration  across  the  out¬ 
let  is  nearly  uniform. 

•  Switching  from  Model_D  to  Model_DC.  This  switch  is 
performed  at  the  moment  when  the  gas  concentration  at 
the  adsorber  wall  becomes  equal  to  the  saturation  concen¬ 
tration  corresponding  to  the  wall  temperature  Csat,  result¬ 
ing  with  the  start  of  condensation.  The  saturation  concen¬ 
tration  is  calculated  using  the  Wagner  equation  (Reid  et 
al.  1987).  The  simulations  show  that  condensation  prac¬ 
tically  starts  at  the  same  time  at  the  whole  surface  of  the 
outer  wall  of  the  adsorber. 

•  Switching  from  Model_DC  to  Model_A.  This  switch  de¬ 
fines  the  end  of  desorption  and  start  of  a  new  adsorption 
part  of  the  cycle.  It  is  performed  when  the  maximal  tem¬ 
perature  in  the  adsorbent  bed  reaches  certain  predefined 
value  Tsw ,  which  should  not  be  exceeded,  and  the  elec¬ 
tricity  is  switched  off.  The  modeled  ESA  cycle  does  not 
assume  a  cooling  step  between  the  end  of  heating  and  start 
of  the  new  adsorption  step  (Sullivan  2003). 

Although  in  principle  different  Comsol  Multiphysics  mod¬ 
els  can  be  run  consecutively,  automatic  switching  from 
one  model  to  another,  by  checking  whether  certain  con¬ 
ditions  are  met,  is  not  possible  in  Comsol  Multiphysics. 
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The  problem  was  solved  by  using  the  coherence  of  Com¬ 
sol  Multiphysics  and  Matlab.  All  three  Comsol  Multi¬ 
physics  models  were  saved  as  Matlab  m-files  and  then 
integrated  into  one  Matlab  model,  defining  the  complete 
cycle.  The  structure  of  this  integral  model  is  shown  in 
Fig.  3. 

In  this  figure,  as  well  as  in  our  program,  the  cycle  is 
started  with  the  desorption  step,  with  an  equilibrium  ini¬ 
tial  state  in  which  the  temperatures  of  both  phases  are  equal 
(and  equal  to  the  temperature  of  the  environment)  and  the 
adsorbent  bed  is  assumed  to  be  saturated,  in  equilibrium 
with  the  inlet  gas  concentration  during  the  adsorption  step. 
This  sequence  was  chosen  because  it  results  with  much 
faster  establishing  of  the  cyclic  steady-state  (in  only  2-3 
cycles),  compared  to  the  case  when  the  process  is  started 
with  the  adsorption  step,  from  an  initially  clean  adsorbent 
bed. 

In  order  to  enable  integration  of  the  models,  the  same 
application  modes  have  to  be  defined  in  all  three  Comsol 
Multiphysics  models,  in  exactly  the  same  sequence.  Also, 
the  model  geometry  and  the  mesh  have  to  be  identical  in  all 
three  models.  Actually,  in  our  integral  model  the  model  seg¬ 
ments  corresponding  to  geometry  definition  and  mesh  ini¬ 
tialization  are  defined  only  once  (Fig.  3). 


5  Solution  of  the  models  and  some  simulation  results 

5.1  Used  mesh,  solvers  and  parameters 

Comsol  Multiphysics  software  is  based  on  the  finite  element 
method  for  solving  systems  of  coupled  partial  differential 
equations  (PDEs)  (Comsol  AB  2005).  The  first  step  in  the 
solution  procedure  is  generation  of  the  mesh  for  the  finite 
element  method.  In  our  simulations  we  used  a  normal  den¬ 
sity  mesh,  with  568  mesh  elements,  which  ensured  a  good 
compromise  between  high  accuracy  of  the  solution,  without 
serious  convergence  problems  on  one,  and  moderate  com¬ 
putational  times,  on  the  other  hand. 

The  finite  element  method  is  actually  a  method  of  dis¬ 
cretization  which  approximates  the  PDEs  by  large  sets  of 
ordinary  differential  equations  (ODEs).  The  software  of¬ 
fers  a  choice  of  different  solvers  for  solving  the  resulting 
sets  of  ODEs.  In  our  simulations  we  used  the  direct  UMF- 
PACK  solver  (Comsol  AB  2005),  which  showed  good  re¬ 
sults. 

The  models  of  the  presented  system  were  solved  and  used 
for  simulation  of  separate  stages  and  of  the  complete  ESA 
cycles,  for  the  following  adsorber  dimensions: 

•  Length  of  the  cartridge:  30  cm 

•  Radius  of  the  central  tube  of  the  cartridge:  0.95  cm 

•  Thickness  of  the  adsorbent  bed:  0.55  cm 

•  Radius  of  the  adsorber  vessel:  3.5  cm. 
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Fig.  3  A  block  scheme  of  the  model  of  the  complete  ESA  cycle  and  the  corresponding  Matlab  program 
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The  physical  parameters  in  our  models  correspond  to  the 
following  adsorption  system  (Sullivan  2003): 

•  Adsorbent:  American  Kynol  ACC-5092-20  (activated 
carbon  fiber  cloth  material) 

•  Adsorbate:  methyl  ethyl  ketone  (MEK) 

•  Carrier  gas:  nitrogen. 

The  simulations  were  performed  for  different  switch  tem¬ 
peratures,  breakthrough  concentrations,  supply  voltages 
and  gas  flow-rates  during  the  adsorption  and  desorption 
half-cycles.  The  inlet  concentration  during  adsorption  was 
equal  to  the  initial  gas  concentration  (in  all  simulations 
0.001  mol/mol).  The  initial,  inlet  gas  and  ambient  tempera¬ 
tures  were  all  equal  (293.15  K). 

The  complete  list  of  the  model  parameters,  with  their  de¬ 
finitions  and  the  numerical  values  used  in  the  simulations  is 
given  in  Appendix  2. 

The  software  offers  a  variety  of  graphical  representation 
of  the  data.  We  will  give  a  small  portion  of  the  simulation 
results,  as  illustration. 

5.2  Simulation  results  for  separate  stages 

The  developed  Comsol  Multiphysics  models  were  first  used 
for  simulation  of  adsorption  on  a  previously  clean  adsorbent 
and  desorption  from  a  previously  saturated  adsorbent  bed. 
As  a  result,  the  gas  velocities,  pressures,  concentrations  and 
temperatures  in  all  parts  of  the  adsorber  are  simulated.  Here 
are  some  of  these  results. 

In  Fig.  4,  the  gas  velocities  in  all  three  subdomains  are 
presented  in  vector  form.  The  velocity  profiles  in  the  system 
are  established  very  fast  (during  the  first  2-3  seconds).  After 


that,  the  velocities  change  somewhat,  owing  to  the  change 
of  the  gas  temperature  and  concentration,  but  the  changes  of 
the  velocity  profiles  are  not  significant.  The  results  shown 
in  Fig.  4  correspond  to  adsorption  and  t  =  7000  s  when  the 
adsorption  is  complete.  Very  similar  picture  of  the  velocity 
profiles  is  obtained  for  desorption.  Figure  4  shows  that  in  the 
central  tube  and  in  the  annular  space  the  gas  flows  mainly 
in  the  axial  direction,  and  through  the  adsorbent  bed  in  the 
radial  direction. 

In  Figs.  5  and  6,  two  samples  of  the  simulation  results, 
representing  the  gas  and  solid  concentrations  and  temper¬ 
atures,  in  the  form  of  surface-plots,  are  given.  The  results 
presented  in  Fig.  5  correspond  to  adsorption,  200  s  af¬ 
ter  the  start  of  adsorption  on  a  previously  clean  adsorbent, 
while  the  results  in  Fig.  6  correspond  to  desorption,  80  s 


Fig.  4  The  velocity  profiles  in 
all  3  sub-domains,  in  vector 
form 
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Fig.  5  Simulation  of  adsorption:  (a)  concentration  of  the  gas  phase;  (b)  concentration  of  the  solid  phase;  (c)  gas  temperature;  (d)  solid  temperature 
( Cin  =  0.001  mol/mol,  Ga  =0.1  mol/s,  t  =  200  s) 
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Fig.  6  Simulation  of  desorption:  (a)  concentration  of  the  gas  phase;  (b)  concentration  of  the  solid  phase;  (c)  gas  temperature;  (d)  solid  temperature 
(Cp  =  0.001  mol/mol,  qp  =  0.00465  mol/g,  U0  =  11  V,  GD  =  0.01  mol/s,  t  =  80  s) 


after  the  start  of  desorption  from  a  previously  saturated  col¬ 
umn. 

By  inspection  of  Figs.  5  and  6  it  can  be  observed  that 
the  concentrations  and  temperatures  change  both  in  the  ax¬ 
ial  and  in  the  radial  direction  of  the  adsorber,  both  in  the 
adsorbent  bed  and  in  the  tubes.  Accordingly,  using  simpler, 
one-dimensional  model,  which  would  neglect  the  changes  in 
the  axial  direction  and  assuming  perfect  mixing  of  the  gas  in 
the  tubes  (Petkovska  et  al.  2005)  would  not  be  justified. 

5.3  Simulation  results  of  the  ESA  cycles 

Simulation  of  the  complete  ESA  process  was  performed 
by  running  the  integral  Matlab  programme,  obtained  by 
integration  the  models  of  the  three  separate  steps,  as  ex¬ 
plained  in  Sect.  4.  A  sample  of  the  simulation  results  of 
the  first  three  ESA  cycles,  corresponding  to  switch  tempera¬ 
ture  463.15  K,  breakthrough  concentration  0.0005  mol/mol 
(50%  of  the  feed  concentration  during  adsorption),  supply 
voltage  1 1  V  and  inert  gas  flow-rates  during  adsorption  and 
desorption  0.1  mol/s  and  0.01  mol/s,  respectively,  is  pre¬ 
sented  in  Figs.  7  (the  outlet  concentration),  8  (the  maximal 
adsorbent  temperature)  and  9  (the  condensation  rate).  All 
other  parameters  were  the  same  as  in  the  simulations  of  the 
separate  adsorption  and  desorption  steps. 

It  can  be  observed  that  after  the  second  cycle  a  quasi¬ 
steady  state  is  practically  achieved.  Also,  very  long  adsorp¬ 
tion  and  short  desorption  times  are  observed. 


Fig.  7  The  outlet  concentration  during  ESA  (C/n  =  0.001  mol/mol, 
U0  =  11  V,  Ga  =  0.1  mol/s,  Gd  =  0.01  mol/s,  Tsw  =  463.15  K  and 
C break  =  0.5  C  in) 

6  Analysis  of  the  influence  of  the  process  parameters  as 
a  basis  for  optimization  of  the  ESA  process 

In  order  to  analyze  the  performance  of  the  ESA  system, 
some  performance  criteria  were  first  defined,  as  measures 
of  the  process  quality: 

•  Adsorption  time,  r  a 

•  Desorption  time,  r  d 

•  Average  outlet  concentration  during  adsorption,  (Cout) A 

•  Average  outlet  concentration  during  desorption,  {Cout) d 

•  Separation  factor — defined  as  the  ratio  of  the  average  out¬ 
let  concentrations  during  adsorption  and  desorption,  as 
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Fig.  8  The  maximal  adsorbent  temperature  during  ESA 
(Cin  =  0.001  mol/mol,  U0=\\Y,Fa=OA  mol/s,  FD  =  0.01  mol/s, 
Tsw  =  463.15  K  and  Cbreak  =  0 .5Cin) 


Fig.  9  The  condensation  rate  during  ESA  (C,-n  =  0.001  mol/mol, 
U0  =  11  V,  Fa  =  0.1  mol/s,  FD  =  0.01  mol/s,  Tsw  =  463.15  K  and 

Cbreak  = 

•  Purification  factor — defined  as  the  ratio  of  the  outlet  and 
inlet  concentrations  during  adsorption,  ap 

•  The  total  amount  of  condensed  liquid  per  cycle,  Lcond 

•  Total  recovery,  defined  as  the  ratio  of  the  total  amount  of 
condensed  liquid  and  the  amount  of  adsorbate  entering 
the  system  during  one  cycle,  Lcondl (Cinz aGa) 

•  Used  electric  energy  per  cycle,  Qei 

•  Used  energy  per  mol  of  condensate,  Qei/Lcond. 

The  influence  of  the  following  operational  parameters  was 
analyzed: 

•  The  switch  temperature  Tsw; 

•  The  breakthrough  concentration  Cbreak \ 

•  The  electrical  voltage  during  desorption  Uq. 


The  results  are  shown  in  tabular  form.  The  influence  of  the 
switch  temperature  on  the  defined  performance  criteria  is 
shown  in  Table  1,  the  influence  of  the  breakthrough  concen¬ 
tration  in  Table  2  and  the  influence  of  the  supply  voltage  in 
Table  3.  All  other  parameters  were  identical  to  those  used  in 
Figs.  7-9. 

The  results  given  in  these  tables  show  clear  influence  of 
the  chosen  operational  parameters  on  the  performance  of  the 
ESA  system.  E.g.,  it  could  be  concluded  that  the  increase  of 
the  switch  temperature  is  favorable  from  the  point  of  view  of 
separation  and  purification  of  the  feed  stream,  as  well  as  re¬ 
generation  of  the  adsorbed  vapor,  but  unfavorable  from  the 
point  of  view  of  energy  consumption.  On  the  other  hand, 
the  increase  of  the  breakthrough  concentration  is  unfavor¬ 
able  from  the  point  of  view  of  separation  and  purification 
of  the  feed  stream,  but  favorable  when  regeneration  of  the 
adsorbed  vapor  and  energy  consumption  are  of  main  inter¬ 
est.  The  increase  of  the  supply  voltage  results  with  higher 
energy  efficiency  (lower  energy  consumption  per  mole  of 
condensate). 

The  total  vapor  recovery  is  ~20%  in  almost  all  cases.  It 
increases  with  the  increase  of  the  switch  temperature,  de¬ 
crease  of  the  breakthrough  concentration  and  with  the  in¬ 
crease  of  the  supply  voltage.  This  value  is  somewhat  un¬ 
derestimated,  as  our  model  assumes  that  condensation  takes 
place  only  during  the  desorption  half-cycle.  Nevertheless, 
the  simulations  showed  that  the  condensation  process  was 
not  finished  at  the  end  of  the  desorption  half-cycle,  i.e.,  when 
the  electricity  was  turned  off.  Consequently,  condensation  is 
expected  to  continue  in  the  beginning  of  the  adsorption-half 
cycle,  until  the  gas  concentration  at  the  adsorber  wall  falls 
below  its  saturation  value.  In  order  to  take  this  into  account, 
another  Comsol  Multiphysics  model  would  need  to  be  in¬ 
corporated,  in  which  the  heat  generation  term  ( 8  Qel/dV )  in 
(20)  is  zero  and  the  boundary  conditions  at  the  adsorber  wall 
correspond  to  condensation  (36). 

7  Conclusions 

In  this  paper  we  have  shown  that  Comsol  Multiphysics™,  in 
combination  with  Matlab™,  can  be  used  for  modeling  and 
simulation  of  a  complex  adsorption  cycle.  Both  the  geome¬ 
try  and  the  phenomena  involved  in  the  modeled  ESA  process 
are  very  complex  (3  sub-domains  and  7  application  modes). 
For  the  chosen  mesh,  the  models  were  used  without  se¬ 
rious  convergence  problems,  and  with  acceptable  comput¬ 
ing  times.  The  integrated  model  of  the  complete  ESA  cy¬ 
cle  was  successfully  used  for  simulation  and  investigation  of 
the  main  operational  parameters  on  the  quality  of  the  ESA 
process. 

The  current  model  was  simplified  in  several  ways.  First, 
most  physical  and  transport  parameters  were  considered  as 
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Table  1  Analysis  of  the  influence  of  the  switch  temperature  Tsw 

Criterium 

Switch  temperature  Tsw  [K] 

433.15 

463.15 

493.15 

(Cout) D  [mol/mol] 

0.00853 

0.00879 

0.00895 

(C0ut) A  [mol/mol] 

2.65  x  10“4 

2.00  x  10“4 

1.84  x  10“4 

&S  =  (CoutM  {Cout)  A 

32.189 

43.950 

48.641 

Mp  =  {C0ut)  a/ Cinf  A 

0.265 

0.200 

0.184 

[mol] 

0.0400 

0.0499 

0.0557 

L cond / (C inT \G \)  [%] 

16.90 

20.45 

21.42 

G*/  [KJ] 

13.363 

16.735 

19.593 

Qell  LCond  [KJ / mol] 

334.07 

335.37 

351.76 

td  [s] 

130 

158 

180 

ta  [s] 

2370 

2440 

2600 

Table  2  Analysis  of  the  influence  of  the  breakthrough  concentration  Cjyreak 


Criterium 

Breakthrough  concentration  Cbreak  [mol/mol] 

2.5  x  1CT4  (25%  C,„) 

3.5  x  1CT4  (35%  Cin ) 

5.0  x  1CT4  (50%  Cin) 

(Cout) D  [mol/mol] 

0.00835 

0.00854 

0.00879 

(Cout) A  [mol/mol] 

1.28  x  10"4 

1.58  x  10"4 

2.00  x  1CT4 

&S  =  {Cout)  d/ {Cout)  A 

65.234 

54.051 

43.950 

®p  =  {C0ut)  a/ Cfn,  a 

0.128 

0.158 

0.200 

Lcond  [mol] 

0.0397 

0.0437 

0.0499 

LCond/ (CinT\G\)  [%] 

22.06 

21.69 

20.45 

Gc/  [KJ] 

15.268 

15.717 

16.735 

Qell LCond  [KJ / mol] 

384.58 

359.66 

335.37 

td  [s] 

143 

148 

158 

ta  [s] 

1800 

2050 

2440 

Table  3  Analysis  of  the  influence  of  the  voltage  used  for  desorption,  Uq 

Criterium 

Supply  voltage  Uq  [V] 

10 

11 

12 

(Cout) D  [mol/mol] 

0.00895 

0.00879 

0.00869 

(Cout) A  [mol/mol] 

1.96  x  10“4 

2.00  x  10“4 

2.08  x  10“4 

=  (Cout)  d/ (Cout)  A 

45.663 

43.950 

41.779 

Mp  —  (C0ut)  a/ Cin^A 

0.196 

0.200 

0.208 

Lcond  [mol] 

0.0500 

0.0499 

0.0510 

LCond/ (CiuTaGa)  [%] 

20.32 

20.45 

20.98 

Qel  [KJ] 

19.626 

16.735 

15.433 

Qell LCond  [KJ/mol] 

392.52 

335.37 

302.61 

TD  [S] 

224 

158 

122 

ta  [s] 

2460 

2440 

2430 
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constants,  although,  strictly  speaking  they  change  with  the 
velocity  and/or  composition  and/or  temperature,  and  con¬ 
sequently,  with  time  and  space  coordinates.  These  rela¬ 
tions  can  be  included  in  the  model,  although  that  would  in¬ 
crease  its  complexity  and  introduce  additional  coupling  of 
the  equations.  Second,  the  heat  capacity  and  resistance  of 
the  adsorber  wall  were  neglected,  although,  they  generally 
should  be  taken  into  account.  In  order  to  do  that,  another 
sub-domain  has  to  be  added,  with  one  application  mode 
(Heat  transfer  by  conduction)  defined  in  it.  Also,  the  heat 
transfer  by  radiation  was  not  taken  into  account,  although 
for  the  temperatures  of  the  adsorbent  bed  during  desorption, 
it  should  probably  not  be  neglected. 

In  the  next  generation  of  multiphysics  models  of  the  ESA 
system  we  are  planning  to  gradually  upgrade  the  complexity 
by  taking  into  account  all  these  aspects. 

Multiphysics  modeling  of  more  complex  adsorbers  with 
two  and  four  cartridges  is  also  in  progress. 
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No.  142014G. 


Appendix  1:  The  model  equations 

8.1  Momentum  balances  and  continuity  equations  for  the 
gas  phase 

For  the  central  tube  and  the  annular  space  between  the  car¬ 
tridge  and  the  adsorber  wall  ( sub-domains  S-l  and  S-3): 


3  u 

rp'i7  +  v 


~  du  ( 

du 

dvY 

_2rM9r  ~r^\ 

Kdz 

+  ^)_ 

3  u 


3  u 


dp 


dv 

rpo - b  V 

8  3 1 


dr  dz 
dv  3  u 


-r'‘l57+fe 


2  r/r 


dv 

dz 


(7) 


dv 


dv 


+  i;— )  +  — ), 


dp 


dz )  dz 


(8) 


3  u  dv\ 


dr  dz  7 


dpg  ,  dpg 


dr 


dz 


:  0. 


(9) 


For  the  adsorbent  bed  ( sub-domain  S-2): 


du 

rp„ - h  V 

Hg  3 1 


du  {  ( du  dv 


M  .  dp\  u 
—u  +  —  I  +  2  ix  — 
k  dr  J  r 


(10) 


dv 

«v>Sj7  +  v 


dv  3  u\  dv' 


IX  dp 

Tv  +  Yz 


The  gas  density  is  calculated  as: 

=  P 
Pg  RgTg(\+Cy 

8.2  Heat  balances  for  the  gas  phase 


(11) 

(12) 

(13) 


For  the  central  tube  and  the  annular  space  between  the  car¬ 
tridge  and  the  adsorber  wall  ( S-l  and  S-3 ): 

r^[Pg(cPg  +  cpvC)Tg ]  +  v(-rZ#^  - 


3C 


3C 


=  -rpgCpyTglu—  +  v—  —  rpg{cPg  +  cpvC)u  — 


dT„ 


dr  dz 
3  Tg 

-rPg(cpg  +  cpvC)v  —  . 

dz 

For  the  adsorbent  bed  (S-2): 

3  r  n 

V  ~dt  ^ 8 ^ Cpg  Cpv ^ Tg  J 


+  VI  —rD"y 


hg  dTg 
7  dr 


rD 


hg  dTg 


t,z,b 


dz 


(  x  (  dC  dC 

'■  r  I  hha(Ts  —  Tg)  —  pgCpVTg  I  u v 

3  Tg 

—  r pg(c pg  +  cpvC)u 

3  Tg 

—  r  pg(c  pg  +  CpVC)  v—^—. 


dr 

(14) 


(15) 


8.3  Adsorbate  balances  for  the  gas  phase 

For  the  central  tube  and  the  annular  space  between  the  car¬ 
tridge  and  the  adsorber  wall  (S-l  and  S-3): 


dC  /  Dmr dC 

r - b  V I  — r - 

3 1  V  Pg 

dC  dC 

=  —  I  ru - b  rv  — 

1  dr  dz 


F>m,z  dC 
Pg  dz 


For  the  adsorbent  bed  (S-2): 


(16) 


dC  (  Dm  r  b  dC  Dm  z  t,  dC 
3 1  \  pg  dr  pg  dz 

kma  (  dC  3 C 

=  r - (C  —  C)  —  I  ru— - b  rv—— 

Pg  \  dr  dz 


(17) 
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8.4  Electric  current  balance  for  resistive  heating  (S-2) 


V 


V  P  dr 


r1-™) 

p  dz  / 


=  0. 


(18) 


Linear  temperature  dependence  of  the  electric  resistivity  of 
the  adsorbent  material,  which  was  obtained  experimentally 
for  ACFC  adsorbent  (Sullivan  2003)  is  used: 


p  =  po(l  +  b(Ts-TR)). 


(19) 


8.5  Heat  balance  for  the  solid  phase  within  the  adsorbent 
bed  (S-2) 


rpb  [(cP,  +  cpIq)Ts]  -  V  ( rD%  ^  +  rD% 

8  Qe,  dq 

=  r~Jv+  rPbAHads-^  ~  rhba(Ts  -  Tg),  (20) 


&Qd 

dV 


1 ffdU\2  fdU^1 

-?((*)  n 


dr  J  ^  \  dz 


1 


p0(l  +  b(Ts  -  Tr))  \\  dr  J  V  dz 


(21) 


8.6  Adsorbate  balance  for  the  solid  phase  within  the 
adsorbent  bed  (S-2) 


dq  . 

PbJ  =kma(C  -  C  ). 

Ot 


(22) 


8.7  Equilibrium  relation 


Fig.  10  Definition  of  the  boundaries 


8.8  Boundary  conditions 


For  better  clarity,  in  Fig.  10  we  show  the  cartridge  geom¬ 
etry  in  which  all  boundaries  between  the  sub-domains,  or 
between  the  adsorber  and  the  surrounding  are  defined  and 
numbered  (from  1  to  13).  The  boundary  conditions  defined 
at  these  boundaries  are: 

Boundary  1 — This  is  the  axis  of  symmetry:  no  momentum, 
mass  or  heat  transfer  through  it. 

Boundary  2: 


Based  on  Dubinin-Radushkevich  equation  (Dubinin  1989), 
a  relation  between  the  equilibrium  adsorbate  concentration 
in  the  gas  phase  C*  and  its  concentration  in  the  solid  phase 
q  and  temperature  Ts  is  obtained: 


RgTs  V  \PaWoJ] 


(23) 


Wagner  equation  (Reid  et  al.  1987)  is  used  to  calculate  the 
adsorbate  saturation  pressure  p°,  needed  in  (23): 


K0 


p°\  _  (VPAx  +  VPBxh5  +  VPcx3  +  VPDxt 


x  =  l-(Ts/Tc). 


(24) 


z  =  0,  r  €  (0,  ri) :  u=  0, 

Tg  =  Tgin,  C  =  Cjfj. 


G 

r\npg ’  (25) 


Boundary  3: 


z  =  H  +  h  i, 
dTk 
dz 


Jig  d 


_  _  n  _  n 

ut,z  M 


r  e  (0,  r\)  :  u  —  0, 

dC 


v  =  0, 


—  =0. 

dz 


(26) 


Boundary  4: 


z  e  (0,  h\), 


■r\  : 


~Dhtf\u  — ^  =  0, 

dz 


dTP 


u  =  0, 

dC  „ 
~Dmr  I  it  ~  —  0. 

dz 


v  =  0, 


(27) 
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Boundary  5: 


Boundary  10: 


r  =  r\,  z£(hi,hi+H):  p\ct  =  p\b. 


u\ct=u\b,  v\ct  =  v\. 


dTn 

~^t,r~dr~  +  Pg(Cpg  +  cpvC)uTg 


~  (  D t§r,b  of  +  Pg(CPg  +CpvC)uTg 


+  hs  i(l  —  £b)(Ts  —  Tg ), 


Dm,r  ~  ~\~  uC 

dr 


—  (  —  Dm,r,b~ - \-uC\  +  km\(l  —  £b)(C*  —  C), 


=hsi(Tg  —  Ts),  J=  0. 

’  or 


Boundary  6 : 


r  =  r2,  z  e  (h,  h  +  H) :  p\b  =  p\t 


U\b  —  V\b  —  V\  as, 


hs  dTn 

~^t,r,b~T~  Pg(CPg  +  cpvC)uTg 


—  ^  Dtfr  Pg(cpg  “I ~  cpvC)uTg 

+  hs  2(1  —  Sb)(Ts  —  Tg), 


dC 

Dm,r,b  ~  h  wC 


—  Dm  r  +  uC  j  +  km2(\  —  &b)(C*  —  C), 

j  =  0,  -Dfrd^  =  hs2(Tg-Ts). 

Boundary  11\ 

z  =  h  +  H,  re  (r2,  7*4) :  1/  =  0,  v  =  0, 


z  —  h\,  re(ri,r2):  w  =  0,  t;  =  0, 


bs  dTg  dC 

-D;*h—  =  0,  -Dm  —  =0, 


t/  =  f/0,  -D%-±=0. 

dz 


hpOtg  at 

-Dt87—=0,  ~Dmz — =0. 

t,z  9z  m.z  9z 

Boundaries  12  and  13 — These  boundaries  correspond  to  the 
outer  wall  of  the  adsorber.  The  boundary  conditions  corre¬ 
sponding  to  this  surface  are  different  for  the  cases  without 
and  with  condensation. 


Boundary  7: 


•  For  Model_A  and  Model_D  (adsorption  and  desorption 
without  condensation): 


z  =  h\  +  H,  r  e  (r\,  r2) :  u  =  0, 

bo  dTg  dC 

~Dt,z,b~T~  =  ~Dm,z,b- ~  =  0, 


U  =  0,  £>*?  —  =0. 

rz  dz 


r  =  r4,  z£(h,h-\-H)  and  re(r2,r4), 

z  G  (0,  /z)  :  u  =  0,  v  =  0, 


3C  „ 
~Dm  r  —  —  0, 
9z 


bo  dTn 

-DZ1^  =  hwt,(Tg-Ta). 


For  Model_DC  (desorption  with  condensation): 


Boundary  8: 


r  =  r4,  zt(h,h-\-H)  and 


z  G  (0,  h\),  r  =  r2:  u  =  0,  v  =  0, 

hs^Tg  dC 

-D%-*-  =  0,  —Dmr—  =  0. 
dz  dz 


r  e  (7*3, 7*4),  z  g  (0,  /*) : 
m  =  0,  v  =  0,  C  =  Csat(Tw), 
hgdTo  dC 

—  ^t,r  =  ~^m,r~^(~2^HCond)  hWg(Tg  —  Ta). 


Boundary  9: 


z  =  0,  r  e  (7*2,  t*3)  :  p\as  =  Pa 


Appendix  2:  Definitions  and  values  of  the  model 
parameters  used  for  simulation 


hodTo  dC 

-Df. — -  =  0,  -Dmz — =0. 

9z  m,z  9z 


All  transport  coefficients  and  physical  parameters  used  in 
the  models  were  considered  as  constants.  The  values  of  the 
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Parameter 


Notation  and  value 


Adsorber  dimensions  (defined  in  Fig.  2) 


Bed  porosity  [5] 

Specific  surface  area  for  heat  and  mass  transfer  [5] 

Ambient  pressure 

Inlet  gas  temperature 

Ambient  temperature 

Wall  temperature 

Bed  hydraulic  permeability  [5] 

Mass  transfer  coefficient  within  the  bed  [1] 

Bed  to  central  tube  mass  transfer  coefficient  [1] 

Bed  to  annular  tube  mass  transfer  coefficient  [1] 

Heat  transfer  coefficient  within  the  bed  [1,2] 

Bed  to  central  tube  heat  transfer  coefficient  [3] 

Bed  to  central  tube  heat  transfer  coefficient  [3] 

Heat  transfer  coefficient  defining  heat  losses  [3] 

Axial  and  radial  dispersion  coefficient  of  the  gas  phase  in  the  central  tube  and  annular  space  [2] 
Axial  and  radial  dispersion  coefficient  of  the  gas  phase  in  the  adsorbent  bed  [2] 

Axial  and  radial  heat  diffusivity  of  the  gas  phase  in  the  central  tube  and  annular  space  [4] 

Axial  and  radial  heat  diffusivity  of  the  gas  phase  in  the  adsorbent  bed  [4] 

Axial  and  radial  heat  diffusivity  of  the  solid  phase  [4] 

Molar  heat  of  adsorption  [5] 

Molar  heat  of  condensation  [5] 

Adsorbent  bed  density  [5] 

Adsorbate  density  [5] 

Adsorbate  molar  mass 
Inert  gas  molar  mass  (nitrogen) 

Specific  heat  capacity  of  liquid  adsorbate  [5] 

Specific  heat  capacity  of  the  inert  gas  [5] 

Heat  capacity  of  the  solid  phase  [5] 

Specific  heat  capacity  of  gaseous  adsorbate  [5] 

Dynamic  viscosity  of  the  inert  gas  [1] 

Electric  resistivity  at  referent  temperature  Tr  =  293  K  [5] 

Temperature  coefficient  of  the  bed  electrical  resistivity  [5] 

Critical  temperature  of  the  inert  gas  (nitrogen)  [5] 

Critical  pressure  of  the  adsorbate  [5] 

Critical  pressure  of  the  inert  gas  (nitrogen)  [5] 

Wagner  constants  for  the  adsorbate  [5] 

Total  volume  of  micropores  (D-R  equation)  [5] 

Adsorption  energy  of  the  adsorbate  (D-R  equation)  [5] 


r\  =  0.95  cm,  r2  =  1 .5  cm 
r3  =  2  cm,  r4  =  3.5  cm 
H  =  30  cm,  h\  =2  cm 
sb  —  0.72 

a  =  13.65  cm2 /cm3 
Fatm  =  101325  Pa 
Tgin  =  293.15  K 
Ta=  293.15  K 
Tw=  293.15  K 
k  =  4.59  x  10“6  cm2 
km  =  0.0003  mol/ (cm2  s) 
km  \  =  0.0003  mol/ (cm2  s) 
km2  =  0.0003  mol/(cm2  s) 
hb  =  0.0055  W / (cm2  s) 
hs i  =  0.001  W / (cm2  s) 
hs 2  =  0.001  W / (cm2  s) 
hwg  =  5  x  10“5  W / (cm2  K) 

Dmz  =  Dmr  =  1  mol/(cms) 

Dm,z,b  =  Dm,z,b  =  8  x  10-5  mol/(cms) 
Df*  =  Dhgr  =  1  W/(K  cm) 

Dhtgz  b  =  Dhgrb  =  0.00014  W/(K cm) 
Dffr  =  Dfr  =  0.0012  W / (cm  K) 

(- A  Hods)  =  66200  J/mol 
(—AHconci)  =  31230  J/mol 
pb  =  0.221  g/cm3 
Pa  =  0.81  g/cm3 
Ma  =  72.107  g/mol 
Mr  =  28.02  g/mol 
cpi  =  157.9  J / (mol  K) 
cpg  =  29.13  J/(molK) 
cPs  =  0.71  J/(gK) 

Cpv  =  100.9  J/(molK) 

/x  =  1.78  x  10-4  g/cm/s 

Po  =  0.202  Q  cm 

b  =  - 2  X  1CT3  1/K 

TcB  =  126.2  K 

PcA  =  4.26  x  106  Pa 

pcB  =  3.39  x  106  Pa 

VP A  =  -7.71476,  VPB  =  1.71061 

VPC  =  -3.6877,  VPD  =  -0.75169 

W0  =  0.748  cm3/g 

E  =  14.43  x  103  J/mol 


Literature  sources:  [1]  Perry  and  Green  (1997),  [2]  Ruthven  (1984),  [3]  Incropera  and  DeWitt  (1996),  [4]  Fuentes  et  al.  (1998),  [5]  Sullivan  (2003) 


transport  coefficients  were  estimated  using  literature  corre¬ 
lations  (defined  by  the  number  in  brackets)  for  some  average 
conditions  during  the  process.  The  values  of  the  physical 


parameters  are  also  taken  from  literature  (given  in  brackets) 
and  correspond  to  some  average  conditions. 
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